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Abstract. Soft Constraint Logic Programming is a natural and flexible 
declarative programming formalism, which allows to model and solve 
real-life problems involving constraints of different types. 
In this paper, after providing a slightly more general and elegant presen- 
tation of the framework, we show how we can apply it to the e-mobility 
problem of coordinating electric vehicles in order to overcome both en- 
ergetic and temporal constraints and so to reduce their running cost. In 
particular, we focus on the journey optimization sub-problem, consider- 
ing sequences of trips from a user's appointment to another one. Solutions 
provide the best alternatives in terms of time and energy consumption, 
including route sequences and possible charging events. 



1 Introduction 

Classical constraint satisfaction problems (CSPs) [12] represent an expressive 
and natural formalism useful to specify different types of real-life problems. A 
CSP can be described as a set of variables associated with a domain of values, and 
a set of constraints. A constraint is a limitation of the possible combinations of 
the values of some variables. So, solving a CSP consists in finding an assignment 
of values to all its variables guaranteeing that all constraints are satisfied. 

Despite their applicability, the main limit suffered by CSPs is the ability 
of just stating if an assignment of certain values to the variables is allowed or 
not. This is indeed not enough to model scenarios where the knowledge is not 
either entirely available or not crisp. In these cases constraints are preferences 
and, when the problem is overconstrained, one would like to find a solution 
that is not so bad, i.e., the best solution according to the levels of preferences. 
For this reason, in |2|3j , the soft CSP framework has been proposed. It extends 
classical constraints by adding to the usual notion of CSP the concept of a 
structure representing the levels of satisfiability or the costs of a constraint. 
Such a structure is represented by a semiring, that is, a set with two operations: 
one (usually denoted by +) is used to generate an ordering over the levels, while 
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the other one (denoted by x ) is used to define how two levels can be combined 
and which level is the result of such a combination. 

Constraint logic programming (CLP) 10 extends logic programming (LP) by 
embedding constraints in it: term equalities is replaced with constraints and the 
basic operation of LP languages, the unification, is replaced by constraint han- 
dling in a constraint system. It therefore inherits the declarative approach of LP, 
according to which the programmer specifies what to compute while disregarding 
how to compute it, by also offering efficient constraint-solving algorithms. 

However, only classical constraints can be handled in the CLP framework. So, 
in [3], it has been extended to also handle soft constraints. This has led to a high- 
level and flexible declarative programming formalism, called Soft CLP (SCLP), 
allowing to easily model and solve real-life problems involving constraints of 
different types. Roughly speaking, SCLP programs are logic programs where 
constraints are represented by predicates which are defined by clauses whose 
body is a value of the semiring modelling the levels of satisfiability or the costs 
of the constraints. The flexibility of the approach is due to the fact that the same 
framework can be used to handle different kinds of soft constraints by simply 
choosing different semirings. It can indeed be used to handle fuzzy, probabilistic, 
prioritized and optimization problems, as well as classical constraints. 

In this paper, before presenting an application of the SCLP framework to the 
e-mobility, we provide a slightly more general and elegant presentation of the 
SCSP framework based on the notion of named semiring, as briefly presented 
in [7]. In particular, soft constraints are elements of the named semiring which, 
besides the additive and multiplicative operations useful to combine the con- 
straints, is equipped by a permutation algebra and a hiding operator allowing 
an explicit handling of names. The permutation algebra indeed allows to char- 
acterize the support of each constraint, that is, the subset of variables on which 
the constraint really depends, while the hiding operator allows us to remove 
variables from the support of a constraint. Then, the SCLP framework can be 
seen as the general case of the SCSP one. 

We propose the SCLP framework as a high-level specification formalism use- 
ful to naturally model and also solve some e-mobility optimization problems. 
With e-mobility indeed new constraints must be considered: electric vehicles 
(EVs) have a limited range and they take long time to charge. So, it is needed to 
guarantee that throughout the user's itinerary the EV never underruns a limit 
energy level and that the user always arrives in time to all his appointments. 

We consider the EV travel optimization problem described in [jj]. In it the 
user has a set of appointments and he makes a series of decisions about the 
sequences of trips from an appointment to another one, for example which route 
to take, if and where to charge the EV and so on. The aim is to find the optimal 
combination of travel choices which minimizes the users cost criteria. In [9] , the 
authors propose a hierarchical presentation of the mobility framework, which 
they exploit to decompose the optimization problem in sub-problems, and in 
particular, they study the journey level one. Here, beside it, we consider the 
optimization problem of the lower level, i.e., the one of the trip level. This last 



problem consists in finding the best trips in terms of travel time and energy 
consumption, while the former one consists in finding the optimal journey, that 
is, the optimal sequence of coupled trips, again in terms of the same criteria, 
guaranteeing that the user reaches each appointment in time and that the state 
of charge (SoC) of the EV never falls below a given threshold. 

The trip level problem substantially coincides with the multicriteria version 
of the shortest path problem modelled in [5] as an SCLP program. So, starting 
from a slightly different specification of this problem, we propose an SCLP pro- 
gram modelling the journey problem. In order to also actually execute both the 
SCLP programs, we propose CIAO Prolog [6], a system supporting constraint 
logic programming. We therefore explicitly implement the soft framework, by 
defining two predicates, the plus and the times ones, which respectively model 
the additive and the multiplicative operations of the semiring. 

This work arises from the research activity in the European project ASCENS 
(Autonomic Service Component ENSembles) aiming at studying formal models, 
languages, and programming tools, for the modelling and the development of 
autonomous, self-aware adaptive systems. E-mobility is right one of the case 
studies of the project and our work represents a first answer to the need of a 
high-level, declarative, executable specification language and of a powerful and 
flexible programming environment where e-mobility problems can be easily and 
naturally modelled and solved. We indeed show as the SCLP framework is able 
to satisfy all these requirements and it can thus be used as a support for rapid 
prototyping and exploratory programming for this kind of problems. 

The paper is organized as follows. Section [5] introduces the SCSP framework 
as an instantiation of the named semiring framework. Section[3]briefly recalls the 
SCLP language and then Section 2] shows how the trip and journey optimization 
problems can be modelled and solved through (S)CLP programs. Finally, Section 
[5] concludes the paper by illustrating some open venues for further works. 

2 Soft Constraints by means of Named Semirings 

This section presents the soft CSP framework based on semiring [213] as an 
instantiation of the more general framework based on named semiring [7] . 

The notion of named semiring is based on the ones of c-semiring (c stands 
for constraint) and permutation algebra. 

2.1 C-Semiring 

Definition 1 (c-semiring). A c-semiring is a tuple (A, +, x,0,l) such that: 

— A is a set and 0, 1 £ A; 

— + : Ax A ^ A is a commutative, associative and idempotent operation, such 
that is its unit element and 1 is its absorbing element; 

— x : A x A — > A is a commutative and associative operation, such that it 
distributes over +, 1 is its unit element, and is its absorbing element. 



Thanks to the idempotence of the + operator, the relation (A,<), defined 
as a < b if a + b = b, is a partial order. Intuitively, a < b means that b is better 
than a or that a implies b. 

It is possible to prove that: (i) the two operations + and x are monotone on 
<; (ii) is its minimum and 1 its maximum; (in) (A, <) is a complete lattice and 
+ is its least upper bouncH. Finally, if x is idempotent then: (iv)+ distribute 
over x; (v) {A, <) is a distributive lattice, and (vi) x is its greatest lower bound. 

2.2 Permutation Algebra 

Here we briefly recall the notion of permutation algebra. We refer the reader to 
[5] for a detailed introduction. 

In the following, we fix a chosen infinite, countable, totally ordered set TV of 
names, which we denote by x, y, z, . . . . 

Definition 2 (Permutations). A name substitution is a function a : M — > 
TV, while a permutation p is a bijective name substitution. The set of all such 
permutations on TV is denoted by 7-* (TV). 

Definition 3 (Kernel). Let p 6 V(Af) be a permutation on TV. The kernel 
of p is the set of the names that are changed by the permutation, formally, 
K(p) = {x eAf\p(x) ^x}. 

A permutation p is finite if its kernel is finite. 

From now on we consider only finite permutations. 

In the following, we introduce the notion of permutation algebra. It consists 
of a pair composed of a carrier set and a description of how the elements of the 
carrier set are transformed by permutations. 

Definition 4 (Permutation Algebras). The permutation signature S p on 
TV is defined as the set of unary operators {p\p £ 'P(TV)} plus the two axioms 
id(a) = a and pi(p2(a)) — (f[p~2{a). 

A permutation algebra A = {A, {pa}) consists of a carrier set A and the set 
of the interpreted operations pa ■ 

Definition 5 (Support). Let A be a permutation algebra and a an element of 
its carrier set. The support of a, supp(a), is the smallest set of names such that, 
given a permutation p, if p{x) = x for all x G supp(a), then Pa(o-) = a. 

Intuitively, supp(a) represents the free names of a: indeed the permutations 
which do not modify them are not influent on a. 

Definition 6 (Finitely supported algebra). A permutation algebra A is 
finitely-supported if each element of its carrier has finite support. 
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2.3 Named c-Semiring. 



A named semiring [7] is a c-semiring plus a finitely-supported permutation alge- 
bra A and a hiding operator (yx.). The permutation algebra allows characterizing 
the finite set of free names of each element c of the named semiring (represented 
by the support of c), while {yx.) applied to c makes the name x local in c. 

Definition 7 (Fusion). A (name) fusion is a total equivalence relation on AT 
with only finitely many non-singular equivalence classes. We denote by x — y 
the fusion with a unique non-singular equivalence class consisting of x and y. 

Definition 8 (Named c-semiring). A named c-semiring C = 
(C, +, x, vx., {pc}, 0, 1) is a tuple where: 

— x = y G C for x, y G Af; 

— (C, +, x,0,l) is a c-semiring; 

— (C, {pc}) is a finite- support permutation algebra; 

— vx. : C —¥ C , for each name x, is a unary operation; 

— for all c,d G C and for all p the following axioms hold: 

(FUSE) x = y X c iff x = y x [y/x]c 

(HIDE) vx.l = 1 vx.vy.c = vy.vx.c vx.{c x d) — c x vx.d if x g" supp(c) 
vx.(c + d) — c + vx.d if x ^ supp(c) vx.c = vy.[y/x]c if y supp(c) 

(PERM) p c = p c l = 1 pc{c x d) = p c c x p c d p c (c + d) = p c c + p c d 
pciyx.c) = vx.(p c c) if x ^ K(p) 

In the (FUSE) axiom, [y/x]c denotes c where y is replaced by x. 
2.4 The Named SCSP Framework 

As briefly shown in [7] , named c-semirings can be suitably instantiated to model 
SCSPs. 

Definition 9 (Constraints). Let S = {A, +, x,0, 1) be a c-semiring, Af a set 
of totally ordered names, and D a finite domain of interpretation for Af. A soft 
constraint is a function (Af — > D) — > A, which associates a value of A to each 
assignment r\ : Af — >• D of the names. 

We define C as the set of all soft constraints over Af, D and A. 

Definition 10. Let S = (A, +, x,0,l) be a c-semiring, Af a set of totally or- 
dered names, and D a finite domain of interpretation for Af. Moreover, let 
C be the set of all soft constraints over Af, D and A. We define the Cscsp 
as the named c-semiring (C, +', <g>, vx., {pc}, 0', 1'), where fusions x — y are 
defined as (x — y)r) = 1 if r)(x) — n(y) and (x — y)r) — otherwise; 
(ci +' c 2 )r\ = cxn + c 2 r\; (ci ® c 2 )i] = cii] x c 2 i]; (vx.c)i] = J2deD( cr l[d/x}); 
(pcc)r/ — erf with w'(x) — i](p(x)); O'rj — and l'r) = 1' for all n. 



The assignment n[d/x] is defined as usual: rj[d/x](y) = d if x = y and 
n[d/x}(y) — rj(y) otherwise. 

Note that, by definition, each constraint c G C involves all the names in Af ', 
but it really depends on the assignment of the names in supp(c) , and intuitively, 
restricting means eliminating a name from the support of the constraint, by 
choosing the best value for it. 

Definition 11 (Soft CSP). A soft constraint satisfaction problem (SCSP) is 
a pair (C, Y) , where C C C is a set of constraints and Y C J\f is a set of names. 

Intuitively, Y represents the set of interface names of the constraint set C . 

Definition 12 (Solution). Let P = (C,Y) be an SCSP. The solution of 
P is the constraint Sol(P) = (vx\) . . .(vx n ){(^C), with {xi,...,X n } = 

Ucec SU PP( C *) \ Y - 

Definition 13 (Best level). Let P = (C,Y) be an SCSP. Then, the 
best level of consistency of P is defined as the constraint blevel(P) — 
(rai)...(ra„)((glC), where {xx,...,x n } = \J Cl ec supp^i)- 

Varying the semiring S, on which the named semiring Cscsp is based, 
several kinds of problems can be represented: we consider the semiring 
Scsp = ({true, false}, V , A, false, true) for classical CSPs; Sfcsp = ({ x \ x G 
[0, 1]}, max, rain, 0, 1) for fuzzy CSPs; and Swcsp = (NU{+oo}, min, +, +oo, 0} 
for optimization CSPs. 

Example 1. Let P be the SCSP with three names TV = {x, y, z}, which can take 
values in D = {red, blue, green}. We assume that only x and y are of interest, 
while the constraints model the fact that all names must take different values. 

We consider the semiring Swcsp, introduced above, and we define three 
constraints, one for each pair of names, c xy , c yz and c zx , associating the worst 
semiring value to the assignments which give the same color to both the names of 
interest for the constraint. In particular, we associate +oo to all the assignments 
giving the same color to both names, if at most one name takes the color red 
then we associate 1, otherwise we associate 2. 

In Horn logic this problem can be expressed as the predicate P{x, y) below: 

P(x,y) : -Q(x,y),Q(y,z),Q(z,x) 

Q{v, w) : — if v = w then +oo else if either v or w are red then 1 else 2. 

The predicate Q(v, w) represents the constraint shown on the left of Fig. (TJ 
where the names v and w represents the only names of the support. Note that in 
this representation we only show the assignment for the names of the support of 
the constraint. Therefore, each table entry actually represents different entries, 
one for each possible color of the names which are not in the support. The con- 
straint representing the problem is then represented by the combination of three 
constraints. Indeed, Q(x, y), Q(y, z) and Q(z, x) respectively represent the three 
constraints c xy , c yz and c zx , obtained by applying to the constraint represented 
by Q(v, w) the permutations mapping the names of the support v and w to the 
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Fig. 1. The constraint represented by Q(v,w) and the solution of the SCSP P. 



names of interest of the three constraints. As said above, the relevant names of 
the problem are just x and y, thus the solution of the problem is obtained by 
combining the three constraints and restricting the scope of the name z. It can 
be expressed as {yz)(c xy ® c yz ® c zx ). The resulting constraint, represented by 
the predicate P(x, y), is shown on the right of Fig. [T] Also in this case, we show 
the assignment of just the two names of the support. Each table entry is the 
minimum value among the ones of the solutions providing a different color for 
z. In this case, the best level of solution is 4, the minimum over all the entries. 

3 Soft Constraint Logic Programming 

This section briefly introduces the soft constraint logic programming (SCLP). 
For a more detailed and complete introduction, we refer the reader to [3]. 

The SCLP framework extends the classical constraint logic programming to 
also handle SCSPs. We can say that an SCLP program over a certain c-semiring 
S is just a CLP program where constraints are defined over S. In the following, 
we fix a semiring S = (A, +, x , 0, 1). 

An SCLP program is hence a set of clauses composed of a head and a body, 
plus a goal. The head of a clause is simply an atom, while the body can be either 
a collection of atoms, or a c-semiring value, or a special symbol □, denoting that 
it is empty. In this two last cases clauses are called facts and define predicates 
representing constraints. When the body is empty, we interpret it as 1, the best 
element of the semiring. Atoms are n-ary predicate symbols followed by a tuple 
of n terms, which can be either a constant or a variable or an n-ary function 
symbol followed by n terms. Ground terms are terms without variables, and 
finally, a goal is a collection of atoms. 

Example 2. As an example, consider the simple SCLP program on the left of 
Fig. [U previously proposed in [3]. We consider the semiring Skcsp and the 
domain D — {a, b, c}. The program is composed of six clauses. The last two are 
facts and the semiring values 2 and 3, associated respectively with the atoms 
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Fig. 2. An example of SCLP program and its fix point semantics. 



t(a) and r(a), mean that they respectively cost 2 and 3 units. The set N of the 
semiring contains all possible costs, and the operations min and + allows us to 
minimize the sum of the costs. We consider as goal the atom :-s(a): later on 
we will show its semantics. 

Three equivalent semantics for the SCLP languages have been defined in 0]: 
the model-theoretic, the fix-point, and the operational one. These semantics are 
conservative extensions of the corresponding ones for logic programming: this 
means that by choosing the c-semiring Scsp we get exactly the LP semantics. 

Actually, we can see the SCLP framework as the general case of the SCSP 
one. As shown in Example [TJ we can indeed express an SCSP program as a set 
of predicates representing the constraints together with a unique clause which 
represents the problem and combines all constraints. In the case of the SCLP, the 
main difference is that the depth of clause nesting is unlimited and the possible 
values D of variables are elements of the Herbrand universe. So, the meaning 
of a predicate P(x, y) assigns a semiring value to all evaluations of variables x 
and y to the Herbrand domain. This is exactly what for example the fix-point 
semantics of the SCLP language does. 

In order to present the fix-point semantics, we need to introduce the notion of 
interpretation and the Tp operator, mapping interpretations into interpretations. 

Definition 14 (Interpretation). An interpretation / consists of a domain D, 
representing the Herbrand universe, together with a function which takes a pred- 
icate and an instantiation of its arguments (that is, a ground atom), and returns 
an element of the semiring: I : {J n (P n —> (D n — ¥ A)), where P n represents the 
set of n-ary predicates and A is the set of the values of the semiring. 

Since interpretations are functions from ground atoms to semiring values, 
we consider programs composed of clauses where the head and the body con- 
tain only ground atoms. So, for example, in the SCLP program of Example [2l 
the clause s(X) :- p(X,Y) is replaced with all its instantiations. In particular, 
for each d G D, we have three clauses s(d) :- p(d,a), s(d) :- p(d,b), and 
s(d) :- p(d,c). 

Definition 15 (Tp Operator). Let P be an SCLP program and LSp the set 

of all its interpretations. Moreover let L be an interpretation and GA a ground 



atom, such that P contains k clauses defining the predicate in GA and the clause 
i is of the shape GA : —B\, . . . , B % n . . Then, we define the operator Tp : ISp — > 

ISp as Tp(I){GA) = X)i=i(Ilj=i Whenever Bj is a semiring value, its 

meaning I(Bj) is fixed in any interpretation I and it is the semiring value itself. 

Definition 16 (Partial Order of Interpretations). Let P be a program and 
ISp the set of all its interpretations. We define the structure (ISp,<), where 
V/1,/2 £ ISp, I\ ^ I2 if It (G A) < I2(GA) for any ground atom GA, where < is 
the order induced by the semiring. 

Note that (ISp, is a complete partial order and its gib coincides with the 
gib operation in the lattice A (extended to interpretations) . Since the function Tp 
is monotone and continuous over this complete partial order, then Tp has a least 
fix-point lfp(Tp) = glb({I\Tp(I) -< I}) and it can be obtained by computing 
Tp j" lo, i.e., by applying Tp to the bottom of the partial order of interpretations, 
and then repeatedly applying it until a fix-point. 

Example 3. Consider again the SCLP program in Fig. [2] In the definition of the 
T P operator we have to consider the additive and multiplicative operations of 
the semiring Skcsp, that is, the min and + operations. As in 0], we start the 
computing of the semantics from the bottom of the partial order of interpreta- 
tions, Iq , which maps each semiring element into itself and each ground atom 
into +00. The table on the right of Fig. [2] shows the value associated by the 
interpretations to the most interesting ground atoms. The interpretation £4 rep- 
resents the fix-point of Tp. As an example, we show the computation for one of 
the most interesting case, the ground atom s(a), which also corresponds to our 
goal. We said that the clause s(X) :- p(X,Y) is considered equivalent to all its 
instantiations. Therefore, £4(3(0)) = min{I^(p(a,a)),I^(p(a,b)),I^(p(a,c))} = 
min{+oo, 2, 3} = 2. 

4 The Electric Vehicle Travel Optimization Problem 

This section presents the EV travel optimization problem, introduced in [9j , and 
shows how it can be naturally modelled and solved in the SCLP framework. 

General description of the problem. A user has a set of appointments, each 
of them is in a location and has a starting time and a duration. The user makes 
a series of decisions regarding the sequences of trips from an appointment to 
another one. For example, he decides which route he wants to follow, where to 
park and if and how to charge the EV at the appointment location. 

All possible combinations of travel choices form the choice set. A travel choice 
is optimal if it minimizes the users cost criteria. 

In particular, finding a single optimal trip consists in finding the best trips 
in terms of travel time and energy consumption. Finding instead the optimal 
journey, that is, the optimal sequence of coupled trips, consists in finding the 
best journeys not only in terms of travel time and energy consumption, but 



also in terms of other important criteria for the user, such as the charging cost, 
the number of charging events, etc. However, in finding the solution it needs to 
guarantee that the user reaches each appointment in time and that the state of 
charge (SoC) of the vehicle never falls below a predefined threshold. 

In [9], the authors propose a hierarchical presentation of the e- mobility 
framework, which they exploit to decompose the optimization problem in sub- 
optimization problems. In particular, they identify four levels of mobility: the 
component level, whose main tasks are the inter- and intra-component coordi- 
nation; the trip level, whose main task is the time and energy optimal routing; 
the journey level, which handles sequences of trips together with charging and 
parking strategies; and the mobility level, which handles mobility services, such 
as, car and ride service. 

Each level represents a different optimization problem and the results of the 
lower level will be inputs of the higher level. However, since, in general, the best 
solution of the lower level could be not optimal for the higher one, the results of 
the lower level could contains several solutions and not only the best one. 

In the following, we consider only the trip level and the journey level opti- 
mization problems. In particular, we first present their formalization and then 
we show how we can model them in the SCLP framework. 

Formalizations of the trip and journey level optimization problems. 

The trip level optimization problem substantially coincides with the multi- 
criteria shortest path problem. The road network is indeed represented by a 
directed graph G := (N, E), where each arc e 6 E from a node p to a node q has 
associated a label (ct,ce), that is, a pair whose elements represent the costs, 
respectively in terms of time and energy consumption, of the arc from p to q. 

So, given the road network G, such as the one on the left of Fig. [3J a source 
node n s and a destination node n^, the problem consists in finding all the best 
paths between n s and rt^ in terms of time and energy consumption. Note that, 
since the costs of the arcs are elements of a partially ordered set, the solution can 
contain several paths, that is, all paths which are not dominated by others, but 
which have different incomparable costs. For example, if we want to know the 
best paths from p to t in the graph of Fig. [31 the solution will contain both the 
paths {p, t} with cost (3,9) and {p, q,t} with cost (4,8). The former is indeed 
better in terms of time while the latter in terms of energy consumption. 

As far as the journey level optimization problem is concerned, we use the 
formalization presented in [5]. Actually, we consider a simpler version of it, which 
avoids to consider car parks and the time that the users would take to go from 
cither the car park or the charging station to the location of the appointment. 
Moreover, we consider only the time and energy consumption as cost criteria 
to be minimized. All these simplifications allows a slender and more readable 
presentation of the SCLP program modelling the problem. 

Let A = {Ai, . . . , A n } be the set of the user's appointments. In order to 
describe the problem, we use different time variables. All of them have the shape 
itX z , where i denotes the appointment, Y G {D,A} (D stands for drive and A 
for appointment), and Z e {S, E} (S stands for start and E for end). 
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Fig. 3. The road network, the user's appointments and the charging stations. 



Each appointments is defined by a location Li, a starting time itg, an end 
time it% and therefore a duration id A . In order to go from an appointment 
Ai to the next one Ai+\, the user leaves with an EV from the location Lj at 
time itg and drive to location Li+%. The user travels along the route alternative 
iR (computed by the trip level problem), which consumes energy and hence 
reduces the SoC. Obviously, the chosen route must allow the user to arrive to 
destination with the SoC of his EV. We assume that the SoC always decreases 
during driving and increases during charging eventfl The user arrives at itg and 
the appointment starts at itg . The user must arrive in time to the appointment, 
so it is required that itg < itg. During the appointment, it is also possible to 
schedule a charging event if the SoC of the EV is not enough to continue the 
journey. We assume to have a set of charging stations. Each of them is simply 
defined by its name CSname, the number of available charging spots SpotsNum, 
and the location L where it is. 

Therefore, given the road network G, a set of appointments, as the ones 
described in the table in the middle of Fig. [31 and a set of charging stations, as 
the ones in the rightmost table of Fig. |3J the problem consists in finding all the 
best journeys through all the appointment locations in terms of time and energy 
consumption. As for the travel optimization problem, also here the solution can 
contain several journeys, that is, all the non-dominated ones. 



SCLP programs for the optimization problems. In the following, we show 
how the SCLP framework can be used as a linguistic support and a high-level 
and flexible programming environment where naturally modeling and solving 
the two optimization problems presented above. 

As far as the trip level optimization problem, we propose a slightly different 
version of the model proposed in [5] for the multi-criteria shortest path problem. 
So, as there, we consider an SCLP program over the c-semiring denoted P H (S) 
which, given a source node n s and a target node rid, allows us to obtain the set 
of the costs of all non-dominated paths from n s to rid- 



4 This is the typical behaviour of EVs, however, as explained in [9], in particular cases 
it might also increase during driving and decreases during charging. 



: -module (paths , _ , _) . 
:-use_module(library(lists)) . 
: -use_module (library (aggregates) ) 

minPair ( [T , E] , [Tl , El] ) : - 
T < Tl, 
E < El. 

times ( [Tl , El] , [T2 , E2] , [T3 , E3] ) : - 
T3 = Tl + T2, 
E3 = El + E2. 

plus([] ,L, []) . 
plus([[P,T,E] IRestL] ,L, 

[[P,T,E] I BestPaths]) :- 

nondominatedC [T,E] , L) , 

plus (RestL ,L ,BestPaths) . 
plus([[P,T,E] IRestL] ,L,BestPaths) 

\+nondominated( [T,E] ,L) , 

plus (RestL ,L ,BestPaths) . 

nondominated ( [T , E] , [] ) . 
nondominated ( [T , E] , [ [P , Tl , El] I L] ) 

\+minPair ( [Tl , El] , [T , E] ) , 

nondominated( [T,E] ,L) . 

Fig. 4. The CIAO program modelling the trip level optimization problem. 



The semiring P H (S) is obtained starting from a semiring S — (A, +, x , 0, 1), 
which in our case is the one modelling the costs associated to each edge, i.e, S = 
(N 2 , min' , +', (oo, oo), (0, 0}), where min' and +' are the min and + operations 
extended to pairs. Indeed, in general we want to minimize the sum of each cost, 
but, since we want to obtain all the non dominated paths, we consider P H (S). 

Given a semiring S, we define P H (S) = (P ff (A),l±l, X*,0, A), where P H (A) 
is the Hoare Power Domain of A, that is, P H (A) = {S C A\x G S,y <s x 
implies j e S}. These sets are isomorphic to those containing just the non- 
dominated values, thus, in the following, we will use this more compact and 
efficient representation, where each element of P H (A) will represent the costs 
of all non-dominated paths from a node to another one. The top element of the 
semiring is the set A (its compact form is {1}, which in our example is {(0, 0)}); 
the bottom element is the empty set; W is the formal union that takes two sets 
and gives their union; x * takes two sets and produces another one obtained by 
multiplying (using the multiplicative operation of the original semiring, in our 
case +') each element of the first set with each element of the second one. 

Note that, in the partial order induced by the additive operation of this 
semiring, a <p»(5) b intuitively means that for each element of a, there exists 
an element of b which dominates it (in the partial order of the original semiring) . 

Following [5 , in order to also really execute the SCLP program, we model 
the problem with a program in CIAO Prolog [BJ, a system supporting CLP, by 
explicitly implementing the soft framework. The program is shown in Fig. 0J 

Here we consider the road network presented in Fig. [3j so we have a set of 
clauses modelling it. In particular, we have a set of facts modelling all the edges 



edge(p,q, [2,4]) . edge(q,t, [2,4]) . 

edge (p , r , [2 , 7] ) . edge (r , s , [3 , 3] ) . 

edge(p,t, [3,9]) . edge(r,q, [1,1]) . 

edge(q,r, [1,1] ) . edge(s,t , [1,1] ) . 
edge(q,s, [4,8]) . 

path(X,Y, [X,Y] ,_, [T,E] ,Lim) :- 
edge(X,Y, [T , E] ) , 
E =< Lim. 

path(X,Y, [XIL] ,V, [T,E] ,Lim) :- 
edge(X,Z, [T1.E1]) , 
nocontainsx(V,Z) , 
path(Z,Y,L, [ZIV] , [T2,E2] ,Lim) , 
times ([Tl, El] , [T2,E2] , [T , E] ) , 
E =< Lim. 

paths (X,Y, Lim, BestPaths) :- 

findall([P,T,E] ,path(X,Y,P, [X] , [T,E] ,Lim) ,ResL) , 
plus(ResL,ResL, BestPaths) . 



Ciao 1.14.2-13646: Moil Aug 15 10:49:59 2011 
?- paths(p,t, 10, BestPaths) . 

BestPaths = [[ [p,t] ,3,9] ,[ [p,q,t] ,2+2,4+4] ] ?. 

no 

?- 

Fig. 5. The output for the query paths(p,t, 10, BestPaths). 

of the graph. Each fact has the shape edge(n s , rid, [ct, c e]), where n s represents 
the source node, rid represents the destination node and the pair [cy, eg] rep- 
resents the costs of the edge in term of time and energy. Note that, differently 
from what would happen in the pure SCLP framework, these facts (represent- 
ing constraints) have the cost in the head of the clauses and not in the body. 
This is needed for implementing the soft framework, and in particular the two 
operations of the semiring. 

Moreover, there are two clauses path describing the structure of paths: the up- 
per one models the base case, where a path is simply an edge, while the lower one 
represents the recursive case, where a path is an edge plus another path. The head 
of the path clauses has the following shape path(n s , rid, Ln , Lv , [ct,ce], Lira), 
where n s and rid are respectively the source and destination nodes, Ln is the list 
needed to remember, at the end, all the visited nodes of the path in the ordering 
of the visit, Lv is the list of the already visited nodes needed to avoid infinite 
recursion where there are graph loops, [ct,ce] is used to remember the cost of 
the path in terms of time and energy, and finally, Lim represents the maximum 
amount of energy that the EV can consume. It is used to retrieve only the paths 
with a total cost in terms of energy equal to or less than the passed value. 

The times and plus clauses are useful to model the soft framework. In par- 
ticular, the first clause is useful to model the multiplicative operation of the 
semiring allowing us to compose the global costs of the edges together, time 
with time and energy with energy. The plus predicate instead mimics the addi- 
tive operation and it is useful to find the best, i.e. non-dominated, paths among 
all the possible solutions. The plus predicate is indeed used in the body of the 
paths clause, which collects all the paths from a given source node to a given 
destination node and returns the best solutions chosen with the help of the plus 
predicate. So, if we want to know the best paths, in the graph of Fig. [3l from 
p to t with a total cost in terms of energy consumption less than or equal to 
10, we have to perform the CIAO query paths{p,t, 10, BestPaths), where the 
BestPaths variable will be instantiated with the list containing all the non- 
dominated paths. In particular, for each of them, the list will contain the se- 
quence of the nodes in the path and the total cost of the path in terms of time 
and energy. The output of the CIAO program for this query is shown in Fig. [5j 

Now, by using the SCLP program modelling the travel optimization problem, 
we can also show the one modelling the journey level problem. Also in this case, 
as before, we consider the P H (S) semiring and we propose a CIAO program, 
where we also model the soft framework. The CIAO program modelling the 
journey optimization problem is presented in Fig. |6l 



: -module C j ourney , _ , _) . 
: -use_module (paths) . 



journey([X,Y] , [P] , [] , [T,E] ,SoC) :- 

appointment CX , Tx , Dx) , appointment (Y , Ty , Dy) , 
path(X,Y,P, [X] , [T,E] ,SoC) , 
timeSum(Tx,Dx,T,ArrT) , ArrT=<Ty. 



plus([] ,L, []). 

plus([[P,T,E,ChEv] IRestL] ,L, 

[ [P , T , E , ChEv] I BestPaths] ) 



journey([X,Y] , [P] , [[X,ID]] , [T,E] ,SoC) :- 

appointment CX , Tx , Dx) , appointment (Y , Ty , Dy) , 
\+path(X,Y,P, [X] , [T,E] ,SoC) , 
chargingStationCID , Spots , X) , Spots>0 , 
newSoC (SoC , Dx , NewSoC) , 
path(X,Y,P, [X] , [T,E] .NewSoC) , 
timeSum(Tx,Dx,T,ArrT) , ArrT=<Ty. 



nondominated( [P,T,E] ,L) , 
plus(RestL,L, BestPaths) . 
plus C[[P,T,E, ChEv] IRestL] ,L, 



BestPaths) 



\+nondominated ( [P , T , E] , L) , 
plusCRestL.L, BestPaths) . 



nondominated ( [P , T , E] , [] ) . 

nondominatedC [P,T,E] , 

[[Pl,Tl,El,ChEvl] |L]) 
\+minPair ( [Tl , El] , [T , E] ) , 
nondominatedC [P,T,E] , L). 



journey([XI [YIZ]] , [P I LP] ,ChEv, [T,E] ,SoC) :- 
appointment CX , Tx , Dx) , appointment (Y , Ty , Dy) , 
path(X,Y,P, [X] , [T1.E1] ,SoC) , 
timeSum(Tx,Dx,Tl,ArrT) , ArrT=<Ty, 
journey([YIZ] , LP, ChEv, [T2,E2] , (SoC-ED) , 
times([Tl,El] , [T2,E2] , [T , E] ) . 



appointmentCp,7,l) . 
appointment Cr , 11 ,2) . 
appointment Ct , 18, 3) . 



chargingStationCcspl , 7,p) . 
chargingStationCcsrl ,4,r) . 
chargingStationCcsr2,0,r) . 



journey ([XI [YIZ]] , [P I LP] , [[X,ID] I ChEv] , [T,E] ,SoC) :- 
appointment (X , Tx , Dx) , appointment (Y , Ty , Dy) , 
\+path(X,Y,P, [X] , [T1,E1] ,SoC) , 
chargingStationCID , Spots , X) , Spots>0 , 
newSoC(SoC,Dx, NewSoC) , 
path(X,Y,P, [X] , [T1,E1] , NewSoC) , 
timeSum(Tx,Dx,Tl,ArrT) , ArrT=<Ty, 
journey([YIZ] , LP, ChEv, [T2,E2] , (NewSoC-ED) , 
times([Tl,El] , [T2.E2] , [T , E] ) . 



journeys (Places , EV.Best Journeies) : - 

f indalK [P,T,E,ChEv] , journey (Places ,P ,ChEv, [T,E] ,SoC) , ResL) , 
plus(ResL,ResL,BestJourneies) . 

Fig. 6. The CIAO program modelling the journey optimization problem. 



We have a set of facts modelling the user's appointments and the charg- 
ing stations. In particular, for each appointment A il there is a clause 
appointment (L i: itg,id A ), while for each charging station we have a clause 
char gingStation(CS name, SpotsNum, L). 

Moreover, there are four journey clauses describing the structure of jour- 
neys. The upper two represent the base case, while the other two represent the 
recursive case. The first clause models the case where a journey is simply a path 
with a cost in terms of energy less than or equal to the SoC of the EV. The 
second clause models the case where the SoC of the EV is not enough to do any 
path and so a charging event, incrementing the energy level, must be scheduled. 
The third journey clause represents the case where a journey is a path with a 
cost in terms of energy less than or equal to the SoC of the EV, plus another 
journey Finally, the last clause models the recursive case where a charging event 
is needed. In all cases we check that the paths allow the user to arrive in time. 

The head of the journey clauses has the shape 
journey(LL,Lp,LchEv,[CT,CE\ 1 SoC), where Ll is the list of the loca- 
tions of the appointments, Lp is the list needed to remember, at the end, all 
the paths of the journey in the correct ordering, LchEv is the list needed to 



Ciao 1.14.2-13646: Moil Aug 15 10:49:59 2011 
?- journeysC [p,r,t] , 10, BestJourneys) . 
BestJourneys = [ 

[[[p,r] , [r,s,t]] ,2+(3+l),7+(3+l), [[r.csrl]]], 
[[[p,r] , [r,q,t]] , 2+C1+2) ,7+(l+4) , [ [r , csrl] ] ] , 
[[[p.q.r] , [r,s,t]] ,2+l+(3+l) ,4+l+(3+l) , [] ] , 
[[[p.q.r] , [r,q,t]] ,2+l+(l+2) ,4+l+(l+4) , []] 
] ?. 

no 

■?- 

Fig. 7. The output for the query journey s(\p,r,t], 10, BestJourneys). 

remember all the charging events needed to complete the journey, [Ct,Ce] 
represents the cost of the journey in terms of time and energy, and finally, SoC 
represents the current energy level of the EV. 

To make the program as readable as possible, we omit the predicates newSoC 
and timeSum, useful to respectively compute the new energy level of the EV 
after a charging event and the arriving time of the user to an appointment. 

The plus clauses are useful to model the soft framework and they are very 
similar to the ones of the trip level problem. The only difference is that here we 
have to consider the charging events. Moreover, note that we reuse the times 
predicate defined in the CIAO program in Fig. [4] 

The journeys clause collects all the journeys through a set of locations (the 
ones of the user's appointments) and returns the best solutions chosen with the 
help of the plus predicate. So, if we want to know the best journeys, in the graph 
of Fig. [21 through the locations where the user has the appointments, with an 
EV having an energy level equal to 10, we have to perform the CIAO Prolog 
query journeys([p, r, t], 10, BestJourneys), where p, r, t are the locations of the 
appointments and the BestJourneys variable will be instantiated with the list 
containing all the non-dominated journeys. In particular, for each of them, the 
list will contain the sequence of the paths of the journey, the total cost of the 
journey in terms of time and energy, and the list of the charging events, each of 
them described by the name of the charging station and its location. The output 
of the CIAO program for this query is shown in Fig. 

5 Conclusion 

In this paper we proposed the SCLP framework as a high-level declarative, ex- 
ecutable specification notation to model in a natural way some aspects of the 
e-mobility optimization problem [9], consisting in coordinating electric vehicles 
in order to overcome both energetic and temporal constraints. In particular, we 
considered the trip and journey optimization sub-problems, consisting in find- 
ing respectively the energy- and time-optimal route from one destination to 
another one, and the optimal sequence of coupled trips, in terms of the same 
criteria, guaranteeing that the user reaches each appointment in time. For both 
the optimization problems, we provided an SCLP program in CIAO Prolog, by 



explicitly implementing the soft framework, that is, the additive and the multi- 
plicative operations of the chosen semiring. The former is a slight variant of the 
CIAO program proposed in [5J Section 4.4] to specify the multicritcria version 
of the shortest path problem. With respect to the program proposed there, here 
we implemented a different semiring, (also proposed in [5j), i.e., the one based 
on the Hoare Power Domain operator, which allowed us to obtain only the best 
(i.e. non-dominated) routes in terms of time and energy consumption. We thus 
provided an implementation of the two operations of this semiring, by defining 
two predicates modelling them. The SCLP program modelling the journey opti- 
mization problem then uses the trip optimization problem results as inputs. It 
is also based on the same semiring that, in this case, allowed us to find the best 
journeys in terms of the two cost criteria. 

As said above, the soft framework is explicitly implemented into each CIAO 
program: there is for example a different plus predicate in each optimization 
program we have proposed. However, it would be interesting to study a general 
way to embed the soft framework in Ciao Prolog. Trivially, one could provide a 
library offering a more general implementation of the operations of the semiring 
of each type of problem. Most interestingly, one could instead think to provide 
a meta-level implementing more efficiently the soft framework. 

Differently from our solution, which allows us to obtain the set of all the opti- 
mal journeys, in the mathematical model proposed in [5] a form of approximation 
is introduced, by considering an aggregated cost function to be optimized. Their 
goal is indeed to minimize this cost function, which considers different cost crite- 
ria: besides the travel time and the consumed energy, they also take into account 
the charging cost, the number of charging events, etc. In modelling the problem, 
here we introduced a simplification by considering just the two main cost criteria, 
that allows us a slender presentation of the work. However, it is obvious that the 
SCLP programs can be easily modified to also take into account several other 
cost criteria. On the other side, we preferred not to introduce any approximation 
of the solution, by instead returning all optimal journeys considered equivalently 
feasible. However, since the use of partially ordered structures, as in our case, can 
in general lead to a potentially exponential number of undominated solutions, 
sometimes it becomes crucial to keep the number of configurations as low as 
possible through some form of approximation allowing us to adopt a total order. 
In this case, the right solution could be to adopt a function that composes all 
the criteria in a single one and then to choose the best tuple of costs according 
to the total ordering that the function induces [5j Section 6.1]. 

As said above, our aim is mainly to propose the SCLP framework as an ex- 
pressive and natural specification language to model optimization problems. We 
indeed think that not only problems representing an extension of the one treated 
here can be modelled by adapting the solution we presented easily enough, but 
that in general our approach can be followed to model hierarchical optimization 
problems. All the SCLP programs we proposed are effective only when data of 
small size are considered. We are indeed conscious that the proposed encodings 
cannot be used to really solve the problem on practical cases, but on the other 



side, we think that CIAO represents a powerful system programming environ- 
ment allowing us not only to write declarative specifications but also to reason 
about them. 

It is therefore clear that here we do not take care of the performances of the 
proposed programs and that our aim is not to compare the performance with 
existing algorithms solving these problems. We indeed leave as future work the 
study of how to improve the performance of our programs. In [SJ Section 8], 
the authors show some possible solutions that could be used towards this end, 
such as tabling and branch-and-bound techniques (implementable for example 
in ECLiPSe [1]). We however would also like to study how our programs can take 
advantage of the use of dynamic programming techniques based, for example, 
on the perfect relaxation algorithm for CSPs [TTj . 

Finally, from a theoretical point of view, as future work, we plan to propose a 
more general framework based on named semiring, allowing us to give a unifying 
presentation of the SCSP and SCLP frameworks, providing an explicit handling 
of the names. 
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